Cancer genomes tolerate deleterious coding mutations through somatic copy number amplifications of wild-type regions.

Cancers evolve under the accumulation of thousands of somatic mutations and chromosomal aberrations. While most coding mutations are deleterious, almost all protein-coding genes lack detectable signals of negative selection. This raises the question of how tumors tolerate such large amounts of deleterious mutations. Using 8,690 tumor samples from The Cancer Genome Atlas, we demonstrate that copy number amplifications frequently cover haploinsufficient genes in mutation-prone regions. This could increase tolerance towards the deleterious impact of mutations by creating safe copies of wild-type regions and, hence, protecting the genes therein. Our findings demonstrate that these potential buffering events are highly influenced by gene functions, essentiality, and mutation impact and that they occur early during tumor evolution. We show how cancer type-specific mutation landscapes drive copy number alteration patterns across cancer types. Ultimately, our work paves the way for the detection of novel cancer vulnerabilities by revealing genes that fall within amplifications likely selected during evolution to mitigate the effect of mutations.

Location and condition based reconstruction of colon cancer microbiome from human RNA sequencing data.

BACKGROUND: The association between microbes and cancer has been reported repeatedly; however, it is not clear if molecular tumour properties are connected to specific microbial colonisation patterns. This is due mainly to the current technical and analytical strategy limitations to characterise tumour-associated bacteria. METHODS: Here, we propose an approach to detect bacterial signals in human RNA sequencing data and associate them with the clinical and molecular properties of the tumours. The method was tested on public datasets from The Cancer Genome Atlas, and its accuracy was assessed on a new cohort of colorectal cancer patients. RESULTS: Our analysis shows that intratumoural microbiome composition is correlated with survival, anatomic location, microsatellite instability, consensus molecular subtype and immune cell infiltration in colon tumours. In particular, we find Faecalibacterium prausnitzii, Coprococcus comes, Bacteroides spp., Fusobacterium spp. and Clostridium spp. to be strongly associated with tumour properties. CONCLUSIONS: We implemented an approach to concurrently analyse clinical and molecular properties of the tumour as well as the composition of the associated microbiome. Our results may improve patient stratification and pave the path for mechanistic studies on microbiota-tumour crosstalk.

Online bias-aware disease module mining with ROBUST-Web.

SUMMARY: We present ROBUST-Web which implements our recently presented ROBUST disease module mining algorithm in a user-friendly web application. ROBUST-Web features seamless downstream disease module exploration via integrated gene set enrichment analysis, tissue expression annotation, and visualization of drug-protein and disease-gene links. Moreover, ROBUST-Web includes bias-aware edge costs for the underlying Steiner tree model as a new algorithmic feature, which allow to correct for study bias in protein-protein interaction networks and further improves the robustness of the computed modules. AVAILABILITY AND IMPLEMENTATION: Web application: https://robust-web.net. Source code of web application and Python package with new bias-aware edge costs: https://github.com/bionetslab/robust-web, https://github.com/bionetslab/robust_bias_aware.

Genes enriched in A/T-ending codons are co-regulated and conserved across mammals.

Codon usage influences gene expression distinctly depending on the cell context. Yet, the importance of codon bias in the simultaneous turnover of specific groups of protein-coding genes remains to be investigated. Here, we find that genes enriched in A/T-ending codons are expressed more coordinately in general and across tissues and development than those enriched in G/C-ending codons. tRNA abundance measurements indicate that this coordination is linked to the expression changes of tRNA isoacceptors reading A/T-ending codons. Genes with similar codon composition are more likely to be part of the same protein complex, especially for genes with A/T-ending codons. The codon preferences of genes with A/T-ending codons are conserved among mammals and other vertebrates. We suggest that this orchestration contributes to tissue-specific and ontogenetic-specific expression, which can facilitate, for instance, timely protein complex formation.

Using protein-per-mRNA differences among human tissues in codon optimization.

BACKGROUND: Codon usage and nucleotide composition of coding sequences have profound effects on protein expression. However, while it is recognized that different tissues have distinct tRNA profiles and codon usages in their transcriptomes, the effect of tissue-specific codon optimality on protein synthesis remains elusive. RESULTS: We leverage existing state-of-the-art transcriptomics and proteomics datasets from the GTEx project and the Human Protein Atlas to compute the protein-to-mRNA ratios of 36 human tissues. Using this as a proxy of translational efficiency, we build a machine learning model that identifies codons enriched or depleted in specific tissues. We detect two clusters of tissues with an opposite pattern of codon preferences. We then use these identified patterns for the development of CUSTOM, a codon optimizer algorithm which suggests a synonymous codon design in order to optimize protein production in a tissue-specific manner. In human cell-line models, we provide evidence that codon optimization should take into account particularities of the translational machinery of the tissues in which the target proteins are expressed and that our approach can design genes with tissue-optimized expression profiles. CONCLUSIONS: We provide proof-of-concept evidence that codon preferences exist in tissue-specific protein synthesis and demonstrate its application to synthetic gene design. We show that CUSTOM can be of benefit in biological and biotechnological applications, such as in the design of tissue-targeted therapies and vaccines.

Single-read tRNA-seq analysis reveals coordination of tRNA modification and aminoacylation and fragmentation.

Transfer RNA (tRNA) utilizes multiple properties of abundance, modification, and aminoacylation in translational regulation. These properties were typically studied one-by-one; however, recent advance in high throughput tRNA sequencing enables their simultaneous assessment in the same sequencing data. How these properties are coordinated at the transcriptome level is an open question. Here, we develop a single-read tRNA analysis pipeline that takes advantage of the pseudo single-molecule nature of tRNA sequencing in NGS libraries. tRNAs are short enough that a single NGS read can represent one tRNA molecule, and can simultaneously report on the status of multiple modifications, aminoacylation, and fragmentation of each molecule. We find correlations among modification-modification, modification-aminoacylation and modification-fragmentation. We identify interdependencies among one of the most common tRNA modifications, m1A58, as coordinators of tissue-specific gene expression. Our method, SingLe-read Analysis of Crosstalks (SLAC), reveals tRNAome-wide networks of modifications, aminoacylation, and fragmentation. We observe changes of these networks under different stresses, and assign a function for tRNA modification in translational regulation and fragment biogenesis. SLAC leverages the richness of the tRNA-seq data and provides new insights on the coordination of tRNA properties.

GZMKhigh CD8+ T effector memory cells are associated with CD15high neutrophil abundance in non-metastatic colorectal tumors and predict poor clinical outcome.

CD8+ T cells are a major prognostic determinant in solid tumors, including colorectal cancer (CRC). However, understanding how the interplay between different immune cells impacts on clinical outcome is still in its infancy. Here, we describe that the interaction of tumor infiltrating neutrophils expressing high levels of CD15 with CD8+ T effector memory cells (TEM) correlates with tumor progression. Mechanistically, stromal cell-derived factor-1 (CXCL12/SDF-1) promotes the retention of neutrophils within tumors, increasing the crosstalk with CD8+ T cells. As a consequence of the contact-mediated interaction with neutrophils, CD8+ T cells are skewed to produce high levels of GZMK, which in turn decreases E-cadherin on the intestinal epithelium and favors tumor progression. Overall, our results highlight the emergence of GZMKhigh CD8+ TEM in non-metastatic CRC tumors as a hallmark driven by the interaction with neutrophils, which could implement current patient stratification and be targeted by novel therapeutics.

Regulation of protein complex partners as a compensatory mechanism in aneuploid tumors.

Aneuploidy, a state of chromosome imbalance, is a hallmark of human tumors, but its role in cancer still remains to be fully elucidated. To understand the consequences of whole-chromosome-level aneuploidies on the proteome, we integrated aneuploidy, transcriptomic, and proteomic data from hundreds of The Cancer Genome Atlas/Clinical Proteomic Tumor Analysis Consortium tumor samples. We found a surprisingly large number of expression changes happened on other, non-aneuploid chromosomes. Moreover, we identified an association between those changes and co-complex members of proteins from aneuploid chromosomes. This co-abundance association is tightly regulated for aggregation-prone aneuploid proteins and those involved in a smaller number of complexes. On the other hand, we observed that complexes of the cellular core machinery are under functional selection to maintain their stoichiometric balance in aneuploid tumors. Ultimately, we provide evidence that those compensatory and functional maintenance mechanisms are established through post-translational control, and that the degree of success of a tumor to deal with aneuploidy-induced stoichiometric imbalance impacts the activation of cellular protein degradation programs and patient survival.

DNA methylation variation along the cancer epigenome and the identification of novel epigenetic driver events.

While large-scale studies applying various statistical approaches have identified hundreds of mutated driver genes across various cancer types, the contribution of epigenetic changes to cancer remains more enigmatic. This is partly due to the fact that certain regions of the cancer genome, due to their genomic and epigenomic properties, are more prone to dysregulated DNA methylation than others. Thus, it has been difficult to distinguish which promoter methylation changes are really driving carcinogenesis from those that are mostly just a reflection of their genomic location. By developing a novel method that corrects for epigenetic covariates, we reveal a small, concise set of potential epigenetic driver events. Interestingly, those changes suggest different modes of epigenetic carcinogenesis: first, we observe recurrent inactivation of known cancer genes across tumour types suggesting a higher convergence on common tumour suppressor pathways than previously anticipated. Second, in prostate cancer, a cancer type with few recurrently mutated genes, we demonstrate how the epigenome primes tumours towards higher tolerance of other aberrations.

Translational adaptation of human viruses to the tissues they infect.

Viruses need to hijack the translational machinery of the host cell for a productive infection to happen. However, given the dynamic landscape of tRNA pools among tissues, it is unclear whether different viruses infecting different tissues have adapted their codon usage toward their tropism. Here, we collect the coding sequences of 502 human-infecting viruses and determine that tropism explains changes in codon usage. Using the tRNA abundances across 23 human tissues from The Cancer Genome Atlas (TCGA), we build an in silico model of translational efficiency that validates the correspondence of the viral codon usage with the translational machinery of their tropism. For instance, we detect that severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is specifically adapted to the upper respiratory tract and alveoli. Furthermore, this correspondence is specifically defined in early viral proteins. The observed tissue-specific translational efficiency could be useful for the development of antiviral therapies and vaccines.

Silencing of SRRM4 suppresses microexon inclusion and promotes tumor growth across cancers.

RNA splicing is widely dysregulated in cancer, frequently due to altered expression or activity of splicing factors (SFs). Microexons are extremely small exons (3-27 nucleotides long) that are highly evolutionarily conserved and play critical roles in promoting neuronal differentiation and development. Inclusion of microexons in mRNA transcripts is mediated by the SF Serine/Arginine Repetitive Matrix 4 (SRRM4), whose expression is largely restricted to neural tissues. However, microexons have been largely overlooked in prior analyses of splicing in cancer, as their small size necessitates specialized computational approaches for their detection. Here, we demonstrate that despite having low expression in normal nonneural tissues, SRRM4 is further silenced in tumors, resulting in the suppression of normal microexon inclusion. Remarkably, SRRM4 is the most consistently silenced SF across all tumor types analyzed, implying a general advantage of microexon down-regulation in cancer independent of its tissue of origin. We show that this silencing is favorable for tumor growth, as decreased SRRM4 expression in tumors is correlated with an increase in mitotic gene expression, and up-regulation of SRRM4 in cancer cell lines dose-dependently inhibits proliferation in vitro and in a mouse xenograft model. Further, this proliferation inhibition is accompanied by induction of neural-like expression and splicing patterns in cancer cells, suggesting that SRRM4 expression shifts the cell state away from proliferation and toward differentiation. We therefore conclude that SRRM4 acts as a proliferation brake, and tumors gain a selective advantage by cutting off this brake.

Protein Complex Organization Imposes Constraints on Proteome Dysregulation in Cancer.

Protein assembly is a highly dynamic process and proteins can interact in different ways and stoichiometries within a complex. The importance of maintaining protein stoichiometry for complex function and avoiding aggregation of orphan subunits has been demonstrated. However, how exactly the organization of proteins into complexes constrains differential protein abundance in extreme cellular conditions like cancer, where a lot of protein abundance changes occur, has not been systematically investigated. To study this, we collected proteomic data made available by the Clinical Proteomic Tumor Analysis Consortium (CPTAC) to quantify proteomic changes during carcinogenesis and systematically tested five interaction types in complexes to investigate which of these features impact on protein abundance correlation patterns in cancer. We found that higher than expected fraction of protein complex subunits does not show changes in their abundances compared to those in the normal samples. Furthermore, we found that the way proteins interact in complexes indeed constrains their co-abundance patterns. Our results highlight the role of the interactions between the proteins and the need of cancer cells to deal with aberrant changes in protein abundance.

Mutation bias within oncogene families is related to proliferation-specific codon usage.

It is well known that in cancer gene families some members are more frequently mutated in tumor samples than their family counterparts. A paradigmatic case of this phenomenon is KRAS from the RAS family. Different explanations have been proposed ranging from differential interaction with other proteins to preferential expression or localization. Interestingly, it has been described that despite the high amino acid identity between RAS family members, KRAS employs an intriguing differential codon usage. Here, we found that this phenomenon is not exclusive to the RAS family. Indeed, in the RAS family and other oncogene families with two or three members, the most prevalently mutated gene in tumor samples employs a differential codon usage that is characteristic of genes involved in proliferation. Prompted by these observations, we chose the RAS family to experimentally demonstrate that the translation efficiency of oncogenes that are preferentially mutated in tumor samples is increased in proliferative cells compared to quiescent cells. These results were further validated by assessing the translation efficiency of KRAS in cell lines that differ in their tRNA expression profile. These differences are related to the cell division rate of the studied cells and thus suggest an important role in context-specific oncogene expression regulation. Altogether, our study demonstrates that dynamic translation programs contribute to shaping the expression profiles of oncogenes. Therefore, we propose this codon bias as a regulatory layer to control cell context-specific expression and explain the differential prevalence of mutations in certain members of oncogene families.

Interactome Mapping Provides a Network of Neurodegenerative Disease Proteins and Uncovers Widespread Protein Aggregation in Affected Brains.

Interactome maps are valuable resources to elucidate protein function and disease mechanisms. Here, we report on an interactome map that focuses on neurodegenerative disease (ND), connects ∼5,000 human proteins via ∼30,000 candidate interactions and is generated by systematic yeast two-hybrid interaction screening of ∼500 ND-related proteins and integration of literature interactions. This network reveals interconnectivity across diseases and links many known ND-causing proteins, such as α-synuclein, TDP-43, and ATXN1, to a host of proteins previously unrelated to NDs. It facilitates the identification of interacting proteins that significantly influence mutant TDP-43 and HTT toxicity in transgenic flies, as well as of ARF-GEP100 that controls misfolding and aggregation of multiple ND-causing proteins in experimental model systems. Furthermore, it enables the prediction of ND-specific subnetworks and the identification of proteins, such as ATXN1 and MKL1, that are abnormally aggregated in postmortem brains of Alzheimer's disease patients, suggesting widespread protein aggregation in NDs.

MIPPIE: the mouse integrated protein-protein interaction reference.

Cells operate and react to environmental signals thanks to a complex network of protein-protein interactions (PPIs), the malfunction of which can severely disrupt cellular homeostasis. As a result, mapping and analyzing protein networks are key to advancing our understanding of biological processes and diseases. An invaluable part of these endeavors has been the house mouse (Mus musculus), the mammalian model organism par excellence, which has provided insights into human biology and disorders. The importance of investigating PPI networks in the context of mouse prompted us to develop the Mouse Integrated Protein-Protein Interaction rEference (MIPPIE). MIPPIE inherits a robust infrastructure from HIPPIE, its sister database of human PPIs, allowing for the assembly of reliable networks supported by different evidence sources and high-quality experimental techniques. MIPPIE networks can be further refined with tissue, directionality and effect information through a user-friendly web interface. Moreover, all MIPPIE data and meta-data can be accessed via a REST web service or downloaded as text files, thus facilitating the integration of mouse PPIs into follow-up bioinformatics pipelines.

Translational efficiency across healthy and tumor tissues is proliferation-related.

Different tissues express genes with particular codon usage and anticodon tRNA repertoires. However, the codon-anticodon co-adaptation in humans is not completely understood, nor is its effect on tissue-specific protein levels. Here, we first validated the accuracy of small RNA-seq for tRNA quantification across five human cell lines. We then analyzed the tRNA abundance of more than 8,000 tumor samples from TCGA, together with their paired mRNA-seq and proteomics data, to determine the Supply-to-Demand Adaptation. We thereby elucidate that the dynamic adaptation of the tRNA pool is largely related to the proliferative state across tissues. The distribution of such tRNA pools over the whole cellular translatome affects the subsequent translational efficiency, which functionally determines a condition-specific expression program both in healthy and tumor states. Furthermore, the aberrant translational efficiency of some codons in cancer, exemplified by ProCCA and GlyGGT, is associated with poor patient survival. The regulation of these tRNA profiles is partly explained by the tRNA gene copy numbers and their promoter DNA methylation.

Generation and Interpretation of Context-Specific Human Protein-Protein Interaction Networks with HIPPIE.

High-throughput techniques for the detection of protein-protein interactions (PPIs) have enabled a systems approach for the study of the living cell. However, the increasing amount of protein interaction data, the varying quality of these measurements, and the lack of context information make it difficult to construct meaningful and reliable protein networks.The Human Integrated Protein-Protein Interaction rEference (HIPPIE) is a web tool that integrates and annotates experimentally supported human PPIs from a heterogeneous set of data sources. In HIPPIE, one can query for the interactors of one or more proteins and generate high-quality and context-specific networks. This chapter highlights HIPPIE's most important features and exemplifies its functionality through a proposed use case.